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Abstract. Coherent Quantum Phase Slip consists in the coherent transfer of 
vortices in superfluids. We investigate this phenomenon in two miscible coherently 
coupled components of a spinor Bose gas conflned in a toroidal trap. After 
imprinting different vortex states on each component, we demonstrate that during 
the whole dynamics the system remains in a linear superposition of two current 
states in spite of the non-linearity and can be mapped onto a linear Josephson 
problem. We propose this system as a good candidate for the realization of a 
Mooij-Harmans qubit and remark its feasibility for implementation in current 
experiments with ®^Rb, since we have used values for the physical parameters 
currently available in laboratories. 


PACS numbers: 03.75.Hh, 03.75.Lm, 03.75.Gg, 67.85.-d 


1. Introduction 

The possibility to set up experiments devoted to the test of quantum phenomena has 
developed a singly peerless increasing interest since the discovery of superconductivity. 
At an early stage, this research led to the realization of Superconducting Quantum 
Interference Devices (SQUIDs) [HU]. Three decades later, the achievement of ultracold 
degenerate quantum gases, like Bose-Einstein condensates (BECs) Oil], opened up 
new opportunities to test quantum interference phenomena, and their implementation 
in Atomtronic Quantum Interference Devices (AQUIDs) (5] E] [7| [H [H [10], the 
atomic analogue of SQUIDs. This duality between superfiuid atomic gases and 
superconductors has stood out both systems as good supports for quantum simulation. 
However, the tunability of the interaction and the versatility of atoms in simulating 
both bosonic and fermionic systems, provide a more promising perspective for the 
implementation of AQUIDs in future technological applications. 

Josephson junctions play a key role in the physics of quantum interference devices. 
They are constituted by two quantum systems connected by a weak link, and can be 
classified in two categories (short and long), owing to the different nature of the 
coupling. In a short Josephson junction, the coherent transfer of physical quantities 
occurs through a single point, the Josephson link mi ns. When the junction is long, 
the coupling occurs locally at each point of the connection. In particular, two spin 
components of a condensate coupled by a Raman laser operate as a long Josephson 
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junction, which is referred in the literature as internal Josephson effect m, and 
obeys the Josephson equations with a coupling proportional to the overlap between 
condensates |T4]. These techniques, by selecting an appropriate spatial dependence of 
the coupling, led to the first observation of vortices in BECs [mnE]. 

Anderson m discussed the role of the phase of the order parameter in superfluids, 
which motivated the study of phase slips in superconductivity [18] , liquid Helium US], 
and BECs jH] |T0] . A phase slip event is a sudden change of the phase in 27r due to 
the motion of quantized vortices through a superfluid. This phenomenon is associated 
to dissipation, as pointed out by Langer and Ambegaokar m- Ultracold atoms, 
as superfluids, can also exhibit phase slips, by winding the phase through solitonic 
states EU m ns m US] and they are able to generate quantum superpositions of 
macroscopic flows. The literature includes several proposals to engineer superpositions 
of flow states in ID quantum gases in continuous rings |26l|271|28| and discrete rings 
El [301 lU [32]. 

Coherent Quantum Phase Slip (CQPS) is an effect recently discovered in 
superconducting systems containing loops [33]. It is the dual phenomenon of the 
Josephson effect, which is a coherent transport of particles between two superfluids, 
but, in contrast, CQPS is defined as the coherent transfer of vortices through the 
Josephson link. In CQPS, the stationary states corresponding to flux quanta states 
become coupled, such that one can continuously change the flux quanta of the system. 
The Mooij-Harmans qubit [34l [35], which consists of a superconducting loop with 
a weak nanowire, was predicted to be able to manifest CQPS between two current 
states. The proposal was made a reality in the experiment of Ref. [33], which led to 
the first experimental observation of CQPS. 

In this work, in order to implement CQPS, we propose the realization of an 
atomic analogue of the Mooij-Harmans qubit by means of a spinor condensate with 
two relevant internal degrees of freedom or spin states. The two components are 
coupled by phase (spin exchange) and density (contact interaction), and both occupy 
the same space region, since the interspecies density repulsion is small enough to keep 
the system in the miscible phase. Therefore, this overlap allows the coupling to occur 
locally, point to point, in the whole bulk of the condensate (long Josephson junction). 
With the aim to engineer a qubit we will select vortex states as the basis of an effective 
two-level system that is able to perform qubit operations [36] . The coherent coupling 
transfers vortices between both components in the absence of population imbalance, 
and the non-linear system exhibits Rabi oscillations. We consider mixtures confined in 
ring geometries, where persistent currents are metastable states and phase slips provide 
the mechanism for the system to exchange winding numbers between components. All 
these properties stand out the system as a promising tool for atomtronie circuits, and 
in particular, for the simulation of CQPS. 

The paper is organized as follows. In sectionj^ we present the mean field model we 
have used to study the system of two miscible coherently-coupled condensates. Section 
characterizes the different dynamical regimes that the system exhibits as a function 
of the coherent coupling and interaction. In section]^ we discuss the regime where the 
system shows CQPS and propose an analytical model that accurately reproduces our 
numerical results. Section [5 is devoted to the other dynamical regimes of the phase 
diagram, and finally, section ^summarizes our work and provides future perspectives. 
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To describe the system of two coherently coupled Bose-Einstein condensates in 
the mean field regime, we will use the Gross-Pitaevskii equation (GPE) for the 
wavefunctions and 

= %v|>^ , (1) 

where Hq = —h?/2mV^ + Etrap? 9i = ^'KhSai/m^ with i is the interaction 

strength, is the scattering length, and Ct is the coherent Raman coupling that forces 
the two components to share the chemical potential ji. We will consider multiply- 
connected geometries by using a toroidal trap of radius R, Etrap = l/2muj‘^p{{p —P 
with angular frequency cjp, and aspect ratio A ^ 1. Eor quasi-2D systems the 
interaction strength takes the value = QiVSirX. 

The time-independent solutions T(r, t) = 0(r) exp (—i/it//i) of Eq. § can be 
found by solving: 

+ 4:Tra^^\4>^\‘^4>^ + 4:Tra^i\(pi\‘^4>^ + 

IJ-h = 'Ho<Pi + 4.'Kaii\(j)if(j)i + 4'Ka^i\(l)-ff(j)i + , (2) 

where all the quantities are written in terms of the harmonic oscillator units, by using 
hujp and ^^h/muJp as energy and length units, respectively. In our simulation we also 
consider = a > to ensure miscibility. The total number of particles N is 

fixed, in such a way that J (j)l{r)(j)i{r)dr = though particles of both components 
can exchange their spin by virtue of the coherent Raman coupling. 

Analytical expressions can be obtained within the Thomas-Eermi approach in 
order to study the ground state properties of the mixture in the regime of large 
interactions. In the case without Raman coupling im El EH, the Thomas-Eermi 
wavefunction is the same for both components: = \/{p — Etrap)/47r(a + 

Such expression will be also useful in the presence of Raman coupling, where the 
ground state satisfies = —ipi [40] for > 0, since a phase difference of tt between 
both components minimizes the mean field energy. Therefore, the Raman coupling Cl 
becomes a simple shift of the chemical potential and the Thomas-Eermi wavefunctions 
can be written as 




I Meff Etrap 

47r(a + 


(3) 


where /ieff = iaPCl/2 will play the role of an effective chemical potential for the coupled 
system. 

We will also explore the effect of vortices on stationary states. Their associated 
angular momentum per particle is hq for each component, where q is the winding 
number or charge of the vortex, when it is centred. This is no longer true for off-centred 
vortices. Nevertheless, one can give an expression for the dependence of the angular 
momentum as a function of the position of the off-centred vortex, by following the 
procedure developed in Refs. miiig for the case of a condensate in a harmonic trap 
in the Thomas-Eermi limit, since in this regime, the effect of the vortex on the density 
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profile can be neglected. We have derived the expression for the contribution on the 
angular momentum per particle of a vortex in the density region of a 2D ring: 


JL 

Nhq 




(1-AA-|(3-A) + 1, 


(4) 


where rj = (vi — R)/6 and RP S are the external and internal Thomas-Fermi radius, 
respectively, and S = = [3(a + /2R\^^^ is the half-width of the torus in 

harmonic oscillator units. 


3. Dynamical regimes of two coupled condensates in a ring 

We investigate the transfer of vortices between two coherently coupled Bose-Einstein 
condensates by solving numerically the time dependent GPE 0 - To this aim, we have 
selected typical experimental values of the physical parameters. We first compute the 
ground state of a two-component ^^Rb spinor condensate with intraspecies scattering 
length a = 101.41 ob and = 100.94 a^, where as is the Bohr radius, corresponding 
to the hyperfine states \F = = 0) and |F = = —!)• The system 

is confined in a toroidal trap, with frequency uOp = 27 t x 200 Hz, aspect ratio 
A = 4, and radius R = 7.5 lam. Afterwards, persistent currents are induced in each 
component, with different winding numbers and by imprinting proper phases. 



Figure 1. Phase diagram containing the different dynamical regimes as a 
function of the effective chemical potential p-eff = n ^ HQ/2 and the Raman 
coupling 12, for = 0.9954, radius R = 7.5/rm, evolving from the initial 

state 1^1 = 1,^2 = 0} in 2D-GPE 0- Solid red line draws the boundary 12c 
between the Trapping regime and the regimes where phase slip exists. Above 12c 
there is a continuous transition from a Non-Coherent Quantum Phase Slip regime 
to a Coherent Quantum Phase Slip regime. The inset compares our numerical 
results for 12c in ID systems (solid curves with open symbols) with the analytical 
expression Eq. 0 (dotted and dashed lines) for the energy gap associated to the 
excitation of spin modes at different values of g^^jg. The labelled points A, B 
correspond to particular cases addressed in later sections. 
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i.e. ^ X ex.p{iqiO), and the system, whose state will be described by the pair 
is let to evolve. 

In Fig. we show the dynamical regimes of the system, for the initial state 
1^1 = 1,^2 = 0), obtained by numerical simulations of Eq. 0 . Depending on the 
values of the effective chemical potential and the Raman coupling, the system explores 
three different regimes that will be explained in detail in the following sections. For 
values of the Raman coupling smaller than a critical one ftc (solid lines) we have found 
a vortex trapping regime, where winding number states can not be exchanged between 
components. We have also explored the influence of the ratio g^ijg by performing 
ID numerical calculations of the spinor GPE, which are shown by solid lines and 
open symbols in the inset of Fig. Open circles correspond to our results for 
g^l/g = 0.9954, whereas triangles correspond to g^i/g = 0.9. 

It has been demonstrated that in order to produce spin excitations [43l |44] , which 
are the relevant ones for phase slips in a spinor condensate, it is necessary to overcome 
an energy gap A given by: 

where n is the total density. For the sake of comparison, we have complemented the 
inset of Fig. [^with the curves given by Eq. 0 for the same numerical values of Qc'. 
dotted {g^i/g = 1.0), dashed {g^i/g = 0.9954) and dot-dashed {g^i/g = 0.9) lines. 
The minimal coupling energy hClc necessary to produce phase slip is of the order of 
A. When g^^ = g, phase slips can be produced for arbitrarily small values of the 
coherent coupling. As the ratio g^i/g decreases the energy cost for producing phase 
slips increases. 

Once this gap is overcome, phase slip is possible. The system continuously transits 
from the Non-Coherent Quantum Phase Slip (NCQPS) regime at O > where 
vortex exchange between components can be observed at rates different from O, to 
the Coherent Quantum Phase Slip (CQPS) regime at ^ We identify a process 
as coherent if it evolves without decay and with a well-defined frequency equal to the 
Raman coupling. 

The dynamical phase diagram also depends on the radius of the torus. As the ring 
geometry constitutes a finite system, a zero-point kinetic energy h? /miC is introduced. 
This energy quantum separates winding number states and, as a result, the degenerate 
states |gi,g 2 ) and |g 2 ,^i) are separated by a gap from other winding number states. 
When the radius of the torus increases, the zero-point kinetic energy goes to zero and 
the energy spectrum becomes a continuum. The same occurs when the interaction 
energy is very large, because the energy to produce a vortex is negligible in front 
of the chemical potential. As we will show later, CQPS decays or is even absent in 
these cases. The dynamical phase diagram can also exhibit dramatic changes in the 
immiscible case, where vortex states can split due to phase separation [45]. 


4. Coherent Quantum Phase Slip 


In order to have an analytical insight into the dynamics of the system, one can follow 
the spirit of the two-mode approximation. In the CQPS regime, the condensate 
wavefunction can be written as usi sg szi: 










V’t,92W 

A,g2 (t) 


(6) 
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where (j)qj{r) = (t)qj{p) x with j = 1,2, are eigenvectors of both the angular 

momentum operator Lz^ with eigenvalue hpj^ and the Hamiltonian without Raman 
coupling, with eigenvalue jUj. This ansatz neglects any contribution from other modes 
with charges different from qi and g 2 - As we will see later, our numerical results agree 
with this assumption, since the only eigenvectors that significantly contribute to the 
dynamics are those associated to the winding numbers imprinted initially onto the 
wavefunction. 

After substituting Eq. ^ in the GPE Q one gets two decoupled linear Josephson 
equations for each winding number: 


and 


ih 

ih 

ih 

ih 


dt 

dt 

dt 

d^l,Q2 

dt 


+ 


w, , 


w, , 

~^Vl,q2 


(7) 


(8) 


The straightforward solution of these linear systems has the eigenvalues /ij ±11/2. 
The energy gap between the two levels is H, which is the driving frequency, then the 
solution for the condensate wavefunction is: 




MiO 


cos 




—i sin ^ 


+ 0^2 {P) ^ 


—i{A^iqt/h—q29+^) 


I sm 




— cos ^ 


,(9) 


and the corresponding densities read: 

1^*1^ = ^{|9^'9i±os2 ^Y^+|(/>qJ^sin2 - Ag6» + (^)},(10) 


where Ag = q 2 — Qi is the initial winding number imbalance, Afiq = 112 — Pi is the 
associated chemical potential imbalance and p is an arbitrary phase. 

Erom Eq. § one can get the mean angular momentum imbalance per particle 
ALz = ^\Lz\^{))/hN^ as a function of time: 

ALz = ^ cos(llt). (11) 

This expression predicts that the exchange of vortices oscillates with the coherent 
coupling frequency H. As a consequence, a 7r-pulse exchanges the winding numbers 
between components, and a tt/ 2-pulse will drive each component to a quantum 
superposition of flows with winding numbers gi and g 2 - 


4 . 1 . Phase slip between adjaeent winding numbers 

Eigurej^ shows our numerical results, within the CQPS regime, for the mean angular 
momentum imbalance per particle, obtained by solving the GPE 0 in condensates 
with different number of particles and O = 0.16cjp (which corresponds to 200 Hz). 
The comparison with the analytical prediction given by Eq. © is also shown. As 
can be seen, the agreement is very good. The frequency of the oscillation of ALz 
is precisely the Raman coupling ft that coherently connects both spin components. 
The population imbalance is initially zero, and remains unaltered during the whole 
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Figure 2. Comparison between the mean angular momentum imbalance per 
particle calculated by solving the GPE for N = b x 10^ (green triangles), 
N = 1.5 X 10"^ (re d squ ares) and N = bx 10^ (blue circles) together with the result 
predicted by Eq. ( |ll| > (thick brown line). The Raman coupling is Cl = 200Hz and 
the initial state is^i = 1, ^2 = 0}- 


a) ^}t=o 





c) Qt=7r 



e) Cit=27T 






Figure 3. Evolution of a condensate with N = b x 10"^ atoms as a function of 
time, after imprinting a vortex on the t component (first row). The second row 
corresponds to the 4- component. The value of the Raman coupling is = 200 Hz 
and the length of the square graphs is 30/rm. The white line corresponds to 
density isocontours at 5% of maximum density, whereas colours represent the 
phase. Panels (a-e) display snapshots of the state during a Rabi cycle. This 
number of particles is in the limit of the CQPS regimes, since ALz deviates in 
6% from the analytical model. 
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simulation, thus one can deduce that the spin exchange occurs at pairs even though 
the population of both spin components is not fixed. 

In order to elucidate how the topological structure of the wavefunction changes as 
a function of time, we show in Fig. |^the dynamical evolution of the density and the 
local phase for a condensate with = 5 x 10^ atoms. White lines represent density 
isocontours at 5% of maximum density and colours depict the phase. Initially, at t = 0, 
we have imprinted a persistent current (gi = 1) in the t component while the other 
is at rest (panel (a)). A quarter of period later, Clt = 7r/2, an azimuthal density node 
[J]is formed spontaneously in each component, at opposite positions (panel (b)). The 
vortex that was inducing the rotation in the initial state escapes from the t component 
through the density depletion, while another vortex crosses the corresponding node 
in the I component, transferring vorticity from the t component to the ^ component 
(panel (c)). This is the mechanism followed by the coupled system to produce 27r- 
phase slips. After that, in panel (d) the evolution is reversed, returning the vorticity 
to the t component (panel (e)), and so on. 

CQPS in atomtronie circuits allows the system to effectively operate as a qubit. 
A quantum mechanical system is a good candidate for qubit manipulation if two 
conditions are fulfilled. First, the system must be considered as an effective two- 
state system, and second, at every time, the qubit must be expressed in a quantum 
superposition of both states. In our system these two states are and 1 ^ 25 ^ 1 ), 

and as a result the state of the system can be written as T = (a|gi,g 2 ) P\Q 2 ^Qi) 
at every time. From the analytical model, we know that a = cos {ftt/2)az and 
P = ishi{VLt/2)(7z^ where cFz is a Pauli matrix. It can be mapped to the most 
general expression of a qubit T = cos(^/2)|0) + exp (i(/:?) sin (^/2)|1). The mapping 
is characterized by a periodic evolution with a period of 27r/^2. At half a period, the 
phases of both components are exchanged, and in between, topological defects appear 
in the wavefunction in order to drive the phase slip. The Raman coupling is a 
parameter that can be externally manipulated, and its control allows to simulate a 
tunable single-qubit quantum gate for quantum information processes. In this regime, 
the system displays two characteristic properties: 

• The system behaves as linear despite the non-linearity, since the Rabi frequency 

is the Raman coupling, independently of the interaction. 

• The evolution occurs following quasi-stationary states. 

The two previous remarkable properties can be demonstrated by following this 
protocol: a) Evolve the system with a certain value of O in the CQPS regime, thus 
the angular momentum imbalance will oscillate with frequency Q. b) At an arbitrary 
time ti, switch the Raman coupling to zero, suppressing the exchange of phase between 
components, c) At another arbitrary time ^ 2 , switch on again the Raman coupling 
to its initial value in the process. Figure represents the whole sequence of this 
protocol. At ti the state gets frozen in a quasi-stationary state rotating at a given 
velocity according to the winding number chemical potential imbalance. Then, at t 2 
the evolution resumes, with exactly the same properties of the system at p. The 
curves predicted by the analytical model accurately fit in with the solution of the 
GPE. 

i This objects should not be confused with dark solitons, since they can also appear in the linear 
case and the associated healing length is, in general, much larger than the one of solitary waves. 
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Figure 4. Comparison between the analytical model given by Eq. |TTJ (thick 
brown lines) and the numerical solution of the GPE (open symbols) for the time 
evolution following the protocol described in the text. Blue squares correspond 
to the mean angular momentum per particle of the t component, and red circles 
to the ^ component. The initial state is \qi = 1,^2 = 0), the number of atoms 
is N = 5 X 10^ and the Raman coupling is = 200 Hz, except between time 
tiQ = 14.37/27r and t 2 ^ = 32.32/27r, where Q is switched off. 


4 . 2 . Phase slip between non-adjaeent winding numbers 

All the results presented in Sect. |4.1| are devoted to the case of CQPS between adjacent 
winding number states. However, the analytical model is more generic and also applies 
between non-adjacent winding numbers. In this section we present the performance 
of CQPS in the case where the winding numbers imprinted onto both wavefunctions 
differ in more than one unit. 

We have shown that in order to change the winding number in one unit, a 27r- 
phase slip event involving the formation of an azimuthal density node has to occur. 
Therefore, to drive each component from winding number qi to q 2 (and viceversa), 
multiple number of such nodes {\qi — q 2 \) must appear simultaneously. Multiple 27r- 
phase slip can not occur through a sequence of single 27r-phase slip events, since other 
states with winding number different from qi and q 2 would contribute, as described 
by the ansatz This fact can be seen in Fig. where we compare the numerical 
results of the GPE (white isocontours at 4% of maximum density and colours for the 
phase), with the analytical prediction for the same density isocontour given by Eq. 
( [Tq| ) (black line), assuming that the system is in the Thomas-Eermi limit and Eq. 
applies, for the initial state \qi = 1,^2 = 0) (panel (a)) and \qi = 2^q2 = 0) (panel 
(b)). The agreement is again excellent. 

Fig. i represents which fixes the critical value of the Raman coupling that 
allows phase slip events, as a function of the effective chemical potential, for different 
initial winding number imbalances /S.q = 1,2,3 (red circles, blue squares and green 
triangles, respectively), after solving the ID-GPE with = 0.9954. Ctc increases 

with /ieff, but this increasing is faster for larger initial /S.q. The azimuthal nodes that 
the system has to generate in order to produce phase slips possess more energy, and the 










Coherent Quantum Phase Slip in two-eomponent bosonie Atomtronie Cireuits 


10 



Figure 5. Comparison of the numerical results of the GPE with the analytical 
prediction of Eq. ( |10| >, for the wavefunction of the t component at a quarter of a 
Rabi cycle. The initial state is \qi = l^q 2 = 0} (panel (a)) and \qi = 2,q2 = 0} 
(panel (b)), the Raman coupling is 12 = 200 Hz and the condensate holds 
N = 5 X 10^ atoms. White lines correspond to density isocontours at 4% of 
maximum density and colours to the phase, both of them obtained numerically. 
Black lines correspond to density isocontours at 4% of maximum density predicted 
by the model, assuming the initial density in the Thomas-Fermi limit, given by 
Eq. 1^. 


strength of the coherent coupling has to be larger to overcome higher energy barriers 
associated to smaller characteristic lengths. 

We have studied the dynamics of the system for different initial winding number 
imbalance. Figure shows the density (white isocontours) and the phase (colour) 
of the t component at a quarter of a Rabi cycle, for different values of the initial 
winding number qi (with q 2 = 0) imprinted onto condensates with N = b x 10^ 
atoms. One can observe that density nodes appear equispaced forming an ordered 
pattern. As the initial angular momentum imbalance increases, the characteristic 



Figure 6. Critical coupling Qc as a function of the effective chemical potential 
)Ueff for different values of the winding number imbalance Aq = 1 (red circles), 
Aq = 2 (blue squares) and Aq = 3 (green triangles), with q 2 = 0. The results 
have been obtained by solving numerically the ID-GPE with g^ilg = 0.9954. 
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Figure 7. Density isocontours at 5% of maximum density and phase (colour) of 
the t component at a quarter of a Rabi cycle for a condensate of N = 10^ atoms, 

12 = 200 Hz and different values of the initial angular momentum imbalance, with 
the 1 component first at rest, a) qi = l,h) qi = 2, c) qi = 3, d) qi = 4, e) qi = 5, 
f) qi = 6, g) qi = 8, h) qi = 10 and i) qi = 16. Only the cases of the first row do 
exhibit CQPS. 


length scale associated to the azimuthal density nodes decreases as predicted by Eq. 
( 10 ), and becomes closer to that of solitary waves as dark solitons or solitonic vortices. 
To generate such objects additional winding number modes have to be excited, and 
therefore, the CQPS process decays after few cycles. Then, the ansatz §) is no longer 
valid. That is why the second and the third row of Fig. will not exhibit CQPS 
(they belong to the NCQPS regime), and only the cases of the first row will display 
this phenomenon. Each solitonic vortex that mediates the phase slip in the cases 
falling in the NCQPS regime contributes to the total mean angular momentum of 
each component according to Eq. 


5. Other dynamical regimes 

As was shown in Fig. CQPS can not be found in the whole parameter space of 
the system, given that g 7 ^ Our numerical results point to the fact that CQPS 
exists as far as the ansatz © is valid, and this occurs for high values of the coherent 
coupling in comparison with Below this critical value, we have found a dynamical 
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regime where there is no vortex exchange between spin components. This fact is due to 
the presence of an energy barrier between winding number states and 1 ^ 25 ^ 1 ), 

that prevents phase slips. Although such a barrier is of nonlinear nature, contrary 
to scalar condensates, it is not due to the existence of solitonic states on the total 
density of the system [25]. As mentioned before, the excitation of spin modes has 
been demonstrated to play a key role in spinor condensates coupled by density, and 
seems to be also relevant in this case. Such excitations are separated from the winding 
number states by the energy gap given in Eq. ([^, which must be overcome in order to 
produce phase slips. If this critical energy can not be transferred between components 
by the coherent coupling, vortices will be trapped, and the mean angular momentum 
per component will not oscillate around the value (gi P q^j2. 



Figure 8. Mean angular momentum per particle of the 1 (solid red line with 
circles) and 1 (dashed blue line with squares) component, after solving the 2D- 
GPE. In the inset, the azimuthal density ani = (a/R) f \(f){p,6)\‘^pdp of the t 
(solid red line), 1 (dashed blue line) component, and the sum of both (dot- 
dashed green line), at fit = 3 (indicated by the dotted vertical line). At the 
right side of the plot, phase pattern of the 1 (top) and 1 (bottom) component 
is represented by colours, and the isocontours at 5% of maximum density by 
the white lines. The effective chemical potential is p-eff = tiWp and the Raman 
coupling is = 2 X 

Fig.|8] shows a typical case representative of the trapping regime. It corresponds 
to the point A indicated in Fig for the coupling = 2 x As can be 

seen, the mean angular momentum of each component oscillates near the initial 
value, and the corresponding densities (shown in the inset after integration along 
the transverse section of the torus) present variations without nodal points. Although 
the interaction between components translate into currents inside each component 
(see the phase maps on the right of the figure), they are not enough to drive phase 
slips. Finally it is worth to note, that during all the time evolution the total density 
remains approximately constant along the torus. 

When the Raman coupling takes intermediate values, CtPCtci stable CQPS will 
not manifest in the dynamics, and the system enters the NCQPS regime. In this 
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case, the coherent coupling is large enough to produce phase slip events that exchange 
the winding number between spin components. However, the time frequency of these 
events is lower than O. This features are reflected in the case displayed in Fig. 
corresponding to the point B of Fig. Now the spin densities can show nodal points 
leading to phase slips, whereas the total density remains constant. Contrary to the 
CQPS case, the position of such nodal points for both components are not located 
at diametrically opposed positions. As the coherent coupling strength increases the 
frequency for phase slips approaches to O. 



Figure 9. Same as Fig. I^for = 5 X 10 ^ujp at Qt = 5.4. 



Figure 10. Density isocontour at 5% of maximum density and phase (colour) 
of the 1 component of a condensate of N = 5 x atoms and Q = 60Hz after 
2 s of evolution, with the state \qi = 1,^2 = 0} as the initial state. Two solitonic 
vortices and a dark soliton appear in the wavefunction. 


In addition, as the chemical potential increases, the excitation of solitary waves 
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(see panels (ci)-(i) in Fig. are responsible for the damping of the exchange of 
angular momentum between components. As a consequence, the system can deviate 
from the quasi-st at ionary path and explore other regions of the phase space. Different 
topological objects are generated, and then, the long-time dynamics will bring the 
condensate to an out-of-equilibrium quantum gas. Notice that in this regime many 
angular momentum modes are excited and the two-mode approach (§ is no longer 
valid. In Fig. [^we show a characteristic snapshot of the density and the phase of the 
t component of a condensate with A" = 5 x 10^ atoms and Raman coupling = 60 Hz 
after 2 s of evolution from the initial state \qi = l^q 2 = 0). The white line traces the 
density isocontour at 5% of maximum density and the colours display the phase. In 
the figure one can also see that although the initial angular momentum imbalance is 
Aq = I, several kinds of solitary waves (a dark soliton and two solitonic vortices) are 
excited in the condensate in order to try to drive a single 27r-phase slip. One can 
see that the appearance of solitonic vortices do not accomplish with the azimuthal 
dependence of the density predicted by Eq. (10), and as a consequence, it is not 
compatible with CQPS. 


6. Summary and conclusions 

In the present work we have proposed an atomic analogue of the Mooij-Harmans qubit 
that displays Coherent Quantum Phase Slip. Two-component condensates loaded on 
toroidal atomtronie circuits can display phase slips by virtue of the coherent coupling. 
When a vortex pattern phase is imprinted onto each component with different winding 
number, the system evolves through quasi-stationary states that are a superposition 
of both winding number states. The two components exchange vortices by phase slip 
events modulated by the coupling, in such a way that the mean angular momentum 
imbalance oscillates with the Raman frequency. 

We have identified the different dynamical regimes of the system as a function 
of the coherent coupling and the effective chemical potential. In particular, we have 
focused on the dynamical phase corresponding to CQPS, where the system behaves 
effectively as linear despite the non-linearity. For this regime, we have mapped the 
dynamics of the coupled system onto linear Josephson equations by using an ansatz 
composed of two winding number modes per component. The whole dynamics, and 
specifically the results obtained for the mean angular momentum imbalance and 
the density, are very accurately reproduced by our analytical model. This model 
predicts that CQPS needs phase slip events to occur through azimuthal density nodes, 
otherwise coherence would be destroyed. Our numerical results obtained by solving 
the time-dependent GPE confirm these predictions. 

We would like to point out the experimental feasibility of this system, since we 
have used values for the physical parameters currently available in laboratories. ^^Rb 
is a good candidate to perform this qubit, mainly, for the closeness of the scattering 
lengths. Toroidal condensates that have been recently obtained have a diameter of the 
order of 10-20 /im, while ours is 15 fim. The Raman coupling does not present any 
strong limitation for its value although commonly, it ranges from tens of Hz, up to 
few hundreds. Besides, phase imprinting techniques have improved in the last decade 
and the individual manipulation of a single component of a spinor BEC is possible 
nowadays. 

The qubit we have proposed points to multiple possibilities in the field of 
cold atoms. The control of the coherent coupling permits to freeze the system 
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in quasi-stationary persistent current states with non-quantized angular momentum 
characteristic of linear superposition of quantum states. Since both current states 
are entangled, one can manipulate the quantum superposition of both flow states, 
performing as a good quantum computer gate, and offering paths for improvements 
in quantum information processing. The theoretical analysis discussing the role of 
non-linear objects, as relative phase domain walls and dark solitons, are out of the 
scope of the present article and will be addressed in the future. 
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